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ANNUAL STATUS REPORT FOR NASA GRANT NAGW-2013 
1. Introduction 

The project entitled “A High Resolution Liquid Xenon Imaging Telescope for 
0.3-10 MeV Gamma-Ray Astrophysics: Construction and Initial Balloon Flights” 
was initiated with the submission in August 1992 of the proposal CAL-2015 in 
response to NASA Announcement of Opportunity NRA 92-OSSA-8. 

For the first year of the project, covering the period 1 January 1993 through 
31 December 1993 funds in the amount of $250,000 were allocated, under NASA 
Grant NAGW-2013. 

The results achieved with a 3.5 liter LXe-TPC prototype during the first year 
can be summarized as follows: 

• The efficiency of detecting the primary scintillation light for event triggering 
has been measured to be higher than 85% (See Section 2.1). 

• The charge response has been measured to be stable to within 0.1% for a period 
of time of about 30 hours. 

• The electron lifetime has been measured to be in excess of 1.3 ms, corresponding 
to attenuation an length of the order of meters. 

• The energy resolution has been measured to be consistent with previous results 
obtained with small volume chambers. At a drift field of 4 kV/cm the FWHM 
of 662 keV gamma-rays is 7% after corrections for electronic noise, shielding 
inefficiency and rise time effects (See Section 2.2). 

• X — Y gamma-ray imaging has been demonstrated with a non-destructive or- 
thogonal wires readout (See Section 2.3). 

• Monte Carlo Simulation Results: detection efficiency, expected background 
count rate at balloon altitude; background reduction algorithms; telescope s 
response to pointlike and diffuse sources; polarization sensitivity calculations 
(See Section 2.4). 

• Design of a 10 liter LXe-TPC laboratory prototype and associated 100 channels 
of readout electronics was completed (See Section 2.5). 

• Construction of a gas purification/recovery system for 10 liter or larger LXe- 
TPC prototype was completed and tested. 

The results obtained in these areas are discussed in details in the following 
sections. They represent a significant progress towards the goal of construction and 
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test of the 10 liter laboratory prototype, with about half the sensitive area and 
depth of the proposed space flight instrument. The work has also addressed and 
answered the questions of technical concern raised by the proposal review panel as 
specified in the proposal CAL-2015 evaluation sheet. 

Three articles published or accepted for publication in refereed journals are 
attached as Appendix 1-3. Three more articles, on which much work was devoted in 
1993, are in the final stage of preparation: the first deals with studies of the different 
background components expected in the LXe-TPC/Coded Mask Telescope in the 
near space environment and their reduction based on imaging (to be submitted to 
Nucl. Instrum. Sc Methods in January 1994); the second deals with the scientific 
use of a LXe-TPC/ Coded mask gamma- ray telescope, with special emphasis on the 
determination of the 26 A1 celestial source distribution (to be submitted to Ap.J. in 
February 1994); the third paper summarizes the results obtained with the 3.5 liter 
LXe chamber triggered by the scintillation light (to be submitted to Nucl. Instrum. 
Sc Methods in March 1994). 

2.1. Triggering on Primary Scintillation Light in LXe 

While the induction signals detected by orthogonal X and Y wire planes give 
two dimensional image of an interaction point inside the detector, triggering with 
the primary scintillation light gives the third coordinate of the interaction point, 
using the known drift velocity. 

The scintillation light trigger efficiency has been studied in great detail with 
our 3.5 LXe liter detector, shown in Fig. 1. It is a gridded ionization chamber with 
a drift gap of 5 cm; the electrode sensitive area is 8 x 8 cm 2 . On the side of the 
chamber body there are two quartz windows with about 70% transmission at the 
wavelength of 175 nm. Two UV sensitive photomultipliers can be fixed on the outer 
cryogenic vessel with vacuum insulation. The distance between the photocathode 
surface of the PMTs and the center of the detector sensitive volume is about 16 cm. 
Two different PMTs have been tested for comparison: an Hamamatsu R2059, and 
an EMI 9813QB. A collimated 137 Cs (10 mCi) 7-ray source is used to irradiate the 
detector 

Each 7-ray interaction inside the detector produces both ionization and scin- 
tillation signals. The scintillation light is detected by the PMTs and used as a 
trigger, while the ionization charge is drifted under the applied electric field and is 
collected by the anode. Fig. 2 shows two typical triggered events, the first one is 
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single photoabsorption event, and the second one is a photoabsorption followed by 
a Compton interaction. 

The waveforms of the light and charge signals are amplified and digitized by a 
CAMAC based LeCroy 2262 FADC, using the light signal as trigger for the FADC. 
Offline data analysis finds the trigger position and the height of the steps of charge 
collection. The steps are fitted with the function, 
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a - pulse height 
b - step starting time position 
c - rise time of the step 
tf - fall time of electronics. 

The delay between the trigger and the time of start of a charge step corresponds to 
the drift time between cathode and grid. 

Fig. 3 shows one typical fitted event. The correlation between step height and 
drift distance is shown in Fig. 4. The dense part in the scatter plot corresponds to 
the full energy peak in the pulse height spectrum. The fact that the peak position 
stays invariant with increasing drift distance, means that there is no charge loss 
during the charge “blob” drift, indicating very good purity of the liquid xenon. 
With the information of the pulse height dependence on drifting distance, one can 
directly infer the attenuation length of electrons in LXe. 

The pulse height distribution gives the spectrum as shown in Fig. 5a. The 
peak is fitted with a Gaussian and a FWHM resolution of 9.8% is obtained at a 
field of 2 kV/cm, which is similar to or better than the resolution of the spectrum 
directly obtained by a MCA, (see Fig. 5b). This shows that the fitting algorithm 
can precisely find the step height, with good noise filtering. It also shows the power 
of the digitized information which permits to reject “bad events” such as those 
interacting in the grid-anode region. 

The geometry of this detector is not optimized for light detection yet. The 
solid angle subtended by the PMTs to the center of the sensitive volume is only 
about 0.5%. The spectrum of scintillation light at zero field and 1.5 kV/cm field 
with 137 Cs light 662 keV 7-ray source is shown in Fig. 6. It shows clearly the effect 
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of light quenching by the external field which reduces the recombination rate. For 
the field higher than 500 V /cm, the fight spectra axe almost the same. This means 
that the quenching effect has reached a saturated value. It can be seen that many 
events produce a single or double photoelectron. However a triggering efficiency 
better than 85% was already achieved for energy depositions larger than 511 keV 
(see Table 1). In Fig. 7, the lower spectrum is from events triggered by a coincidence 
of two PMTs, while the upper is the spectrum from direct charge signal. From the 
ratio of the total events under the peak, one gets the trigger efficiency of 73.3% for 
two PMT coincidence triggering. Random coincidences, which are less than 0.1% 
of the total counts, were subtracted when estimating the triggering efficiency. 

Table 1. Triggering efficiency for different energy fines 
with single PMT 


511 keV 662 keV 1.274 keV 
Trig. Eff. (%) 86.9 87.2 100 


The triggering efficiency is expected to be improved for the 10 liter LXe-TPC 
prototype due to the better geometry for fight detection. 

2.2. Charge and Energy Resolution Response 

The same 3.5 liter chamber shown in Fig. 1, with a drift distance of 5 cm, was 
used for studies of charge collection and energy resolution dependence on: gamma- 
ray energy, shielding effects, pulse rise time variation, and source collimation. Ex- 
cellent LXe purity was achieved for the 3.5 liter sample, as demonstrated by the 
ionization pulses of cosmic ray muons traversing the 5 cm gap and recorded at a 
filed of 5 V/cm and 2 V/cm (Fig. 8). The drifting electron lifetime is in excess of 
1.3 ms, which is the longest pulse rise time recorded at the lowest field. 

Charge collection stability and energy resolution was measured with a colli- 
mated 137 Cs source, as well as with 22 Na and 207 Bi sources. The relevant results 
obtained with 662 keV gamma-rays are shown in Fig. 9(a-d). Fig. 9a shows the 
charge yield as a function of applied electric field; Fig 9b shows the dependence of 
the noise subtracted energy resolution on the drift field. At 5 kV/cm the FWHM 
resolution at 662 keV is 8%, in good agreement with previous results obtained at 
the same field with small volume LXe chambers. The degradation of resolution 
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due to shielding inefficiency and pulse rise time variation is not negligible in the 
present data. The contribution was estimated by Monte Carlo. When subtracted 
in quadrature, along with the electronic noise contribution, the intrinsic resolution 
becomes on the order or 7% FWHM at 662 keV. Fig. 9c shows the measured res- 
olution dependence on gamma-ray energy. A good fit with a 1 /y/~E function is 
obtained. The data were obtained at a drift field of 2.5 kV/cm. Typical 137 Cs 
and 22 Na spectra obtained at low fields are shown in Fig. 10. Fig 9d shows the 
variation of the collected charge as a function of time. Over a period of about 30 
hours the detector’s response was stable to within 0.1%. This demonstrates that 
there should be no concerns about the lifetime of LXe against buildup of impurities 
due to outgassing. As for radiation damage in the space environment, it is known 
that one of the advantages of liquid rare gases is their radiation hardness. 

2.3. X-Y Gamma-Ray Imaging 

In the few MeV region, Compton interaction dominates and produces events 
of multiple interaction sites in the LXe 7-ray detector. These multiple-site events 
are characterized by the spatially separated point-like charge blobs. These might 
give the ambiguity of event localization and degrade the resolution of imaging. 
Therefore it is necessary to accurately localize these charge blobs and precisely 
measure the energy deposited in each site, in order to reconstruct the event based 
on the Compton kinematics. 

We adopted in our proposal a simple non-destructive readout, schematically 
shown in Fig. 11a, for imaging the point-like charge blobs produced by 7-ray inter- 
actions in LXe. It consists of two orthogonal induction wire planes, separated from 
the drift region by a screening mesh. The sense wires have a spacing that is a mul- 
tiple of the mesh pitch and are aligned with the mesh so that they don’t intercept 
the drifting electrons. The induced signals on the wires provide X — Y information. 
The measured drift time, referred to the scintillation trigger, and the known drift 
velocity, provide the Z- information. Thus 3-D imaging is achieved for each ionizing 
event. The charge of each blob can be obtained from the multiple-steps pulse shape 
recorded by an anode plate placed below the induction wires. 

Such a readout geometry turns out to be very suitable for accurate localization 
of the point-like charge blobs. The amplitude of the induction signal produced 
by a point-like charge depends strongly on its initial lateral position with respect 
to the sense wire cell. This position dependence of the induction signals can be 
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an advantage to localize a point-like charge blob by weighting the signals from 
adjacent wires. This is schematically shown in Fig. lib. Without the weighting, the 
spatial resolution is limited by the wire spacing. By pulse height sharing between 
adjacent wires and subsequent centroid finding, the position of a charge blob is 
determined with an effective resolution < 58/zm(2mm/\/l2) for 2 mm wire spacing. 
The ultimate resolution is therefore determined by the accu- racy of the amplitude 
measurements, i.e. the signal- to- noise ratio. 

We have successfully tested this design concept with the 3.5 liter LXe chamber. 
The shielding mesh is made of nickel, spot-welded on a stainless steel frame, and 
has a spacing of 0.8 mm and a wire width of 50 /jm. The ratio of electric fields on 
the two sides of the mesh was chosen to be 2.0 or more to ensure full transparency 
for the drifting electrons. The sense wires are made of BeCu with a diameter of 50 
mm and have a spacing of 4 mm, which is a multiple of the mesh pitch. The tested 
readout has 10 sense wires for the X-direction and 10 for the F-direction. The 
sense wires are stretched on a ceramic (MACOR) holder, electrically isolated from 
each other and biased with their natural potentials in order not to distort the field 
distribution. The distances between the plates, determined by the ceramic spacers, 
are G x — 3.25 mm, G y = 3.75 mm and G a = 4.0 mm (see Fig. 1). The induction 
signal on each wire is detected, via a decoupling capacitor, with a low noise charge 
sensitive preamplifier with a decay time of 50 ms. The output of each preamplifier 
is connected to a FADC through an additional amplification stage. This additional 
amplification stage has a band-pass filter which greatly enhances the signal-to-noise 
ratio of the induction signals without changing the triangle shape of the signal. 

Fig. lib shows a typical X — F image of a 7-ray event from 207 Bi, recorded 
by the FADC readout at a drift field of 0.5 kV/cm. Due to the point-like charge 
blob produced by the 7-ray, only two adjacent X- wires and two adjacent F- wires 
show a signal. The different pulse heights on these adjacent wires give better lo- 
calization accuracy than that of wire spacing, by weighting the signal amplitudes. 
This is further demonstrated in Fig. 12 by the pulse height distribution of the in- 
duction signals which are induced on a particular X-wire by blobs of fixed amount 
of charge. The blob with fixed amount of charge was selected by setting a nar- 
row window on the energy spectrum and picking up events with single step pulse 
shape. These charge blobs are produced by the 7-rays of fixed energy (662 keV) 
and uniformly interacting within the drift region (uncollimated source). The two 
bumps correspond to two bunches of possible drifting paths for the charge from 
the drifting region, which are focused by the shielding mesh, as shown in Fig. 12. 
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The whole distribution corresponds to the lateral distance away from the monitored 
wire within the 4 mm pitch for X- wires. The maximum and minimum pulse heights 
correspond to the interaction point which is laterally closest or farthest from the 
wire, within a wire cell. The RMS spread of these bumps directly gives the spatial 
resolution for this horizontal axis, which is obviously better than the wire spacing, 
on an event-by-event basis and is approaching the intrinsic limit set by the physical 
extension and diffusion of the charge blobs. The step height of the anode signal 
gives the total deposited energy. 

In this preliminary design and test, we have demonstrated that both 7-ray 
imaging with submillimeter X — Y localization accuracy and calorimetry can be 
obtained simultaneously, using such a readout. Fig. 13 shows a distribution of the 
events when imaging collimated 7-rays from 137 Cs with a beam width of about 4 
mm. 


2.4. Monte Carlo Results 

MonteCarlo calculations continued throughout 1993, with the goal of finding 
the optimum detector and shield configuration for best sensitivity in the energy re- 
gion of 0.3-10 MeV. The capability of the LXe gamma-ray detector to measure both 
the deposited energy and the three spatial coordinates of any interaction point inside 
the sensitive volume, enables one to use the kinematics of the dominant Compton 
scattering process to identify the direction of a multiple-sited 7-ray event and there- 
fore translates directily into a powerful background rejection capability. The event 
reconstruction algorithm, based on Compton kinematics and originally discussed 
in the published NIM paper (Appendix 1), was futher tested and optimized for 
different modes of telescope operation. Algorithm 1 was developed for operation 
with a fine resolution mask with the 7-ray source located within a narrow FOV. 
In this case the initial 7-ray direction is assumed perpendicular to the detector’s 
plane. Algorithm 2 was developed for a wide FOV telescope, optimized for diffuse 
sources observations. Both algorithms are most efficicient in the energy range 0.5-3 
MeV, where Compton cross section is dominant. All events with a single point of 
interaction are rejected as background and no pair production analysis is done to 
recover high energy events. The photopeak effiency of the LXe-TPC/Coded Mask 
Telescope was calculated for the energy range 0.3-10 MeV for 5 cm (laboratory 
prototype) and 10 cm (flight prototype) active depth of LXe. Fig. 14 shows the 
efficiency for 5 cm, before and after the Compton event reconstruction. To recover 
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the efficiency at low energy, more sophisticated analysis based on energy and loca- 
tion can be applied. Similarly high energy events characterized by pair production 
can be recovered by appropriate analysis. The effect of the event reconstruction on 
Compton continuum suppression is shown in Fig. 15. To estimate the background 
rate expected in LXe-TPC/Coded Mask Telescope at balloon altitude (3 g/cm 2 over 
Palestine), we considered the contribution from atmospheric and cosmic 7-rays, as 
well as background originating from the detector itself and the materials around the 
detector. Results from the calculations are shown in Figs. 16-20. The calculations 
were done with either 5 cm or 10 cm thick active Csl shield around the sides and 
bottom of the LXe-TPC. The figures show the total background count rate as well 
as the components considered: aperture flux, shield leakage, mask, spallation and 
elastic neutron scattering. The reduction of the background count rate by the event 
reconstruction algorithms is also shown. For spallation induced background the re- 
duction efficiency was assumed to be the same as that for atmospheric background. 
In reality the reduction factor should be even larger because most of the spallation 
events will be beta-decays. The background rates for different configurations are 
summarized in Table 2. The large contribution at low energies (0. 1-0.3 MeV) will 
be strongly reduced by the scintillation trigger threshold as well as by the plastic 
scintillation counters to be used for charged particles veto in front of the aperture. 
Results indicate that a shielded LXe instrument will not be dead-time limited, even 
with a maximum event rate of about 600/s for the entire 0.1—10 MeV region. Calcu- 
lations with an unshielded detector show that the event rejection technique applied 
to an unshielded detector can provide the same reduction factor as a detector with 
thick shield. Details of the background components studies including neutron acti- 
vation and the reduction techniques based on event imaging are summarized in a 
paper in preparation for NIM. Here we discuss in particular the secondary radiation 
from the tungsten coded mask which was considered of particular concern by the 
proposal review panel. 

Table 2. Total Background Rate of LXe Telescope 
at 3 g/cm 2 Over Palestine, Tx 


Before rejection 


After rejection 5 cm shield 


0.1-10 MeV 
0.3-10 MeV 


10 cm shield 
580 s" 1 
180 s" 1 


5 cm shield 
630 s- 1 
230 s" 1 


Algorithm 1 
70 s- 1 
70s- 1 


Algorithm 2 
50 s- 1 
50 s- 1 
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Gamma-ray background may be produced in the mask from (1) ambient 7-rays, 
(2) charged particles (photons and electrons), and (3) neutrons. The first source, 
ambient photons, can Compton scatter or produce pairs in the mask material. The 
scattered photons and bremsstrahlung photons (from electromagnetic cascade, s for 
example) can then interact in the detector. Using the empirical formula of Owens 
et al. (1985) and the count rate spectrum measured in a balloon flight by the UNH 
group (Dunphy et al. 1989), we find that most of the mask background is due to 
reprocessing of high energy ambient (i.e., atmospheric) photons. 

The second source, 7-rays produced by electrons and protons (> 20 MmeV) 
incident on the mask, can be estimated from measurements and calculations of 
electron and proton fluxes at balloon altitudes (Verma 1967; Daniel and Stephens 
1974; Aitbaev et al. 1981). This can be compared with the atmospheric 7-ray 
flux (Thompson 1974; Schonfelder et al. 1980). This shows that the flux of 7-rays 
produced in the mask by charged particles is at least a factor of 5 smaller than 
that produced by the atmospheric photons, and thus can be neglected. The third 
source, background photons from the interaction of thermal and fast neutrons in the 
atmosphere with a tungsten mask, has been evaluated by Butler et al. (1984) using a 
Monte Carlo code. Using their calculated yields, we find that the photon background 
from this source is about 4 orders of magnitude smaller than the atmospheric 7-ray 
background from the region of sky covered by the mask; so this is also negligible. 

To calculate the contribution of the mask background to the LXe detector 
counting rate, we compare it to the background from the atmosphere and the cosmic 
diffuse gamma-ray flux. Comparing the count rate spectrum from the mask in 
Dunphy et al. (1989) at 1 MeV with the count rate expected from the atmosphere 
plus the cosmic diffuse source (Schonfelder et al. 1980; Costa et al. 1984), we 
conclude that the mask adds a relatively small component (50% or less) to the 
photon background that comes through the aperture of a coded aperture detector. 
The estimated background counting rate spectrum in the LXe detector from the 
mask is shown in the previous Figs. 16-20. 

Using the calculated efficiency and background rate after event reconstruction, 
we have simulated the response of the telescope to various source distributions. To 
illustrate the point source response we show in Fig. 21 the image of the Crab Nebula 
in the range 1-2 MeV (10 hr observation). The source is present at a significance level 
of 10 s. The localization accuracy for this observation is 3 areminutes, assuming a 
characteristic pixel size of 30 areminutes. The performance of the proposed telescope 
as a Compton Polarimeter is discussed in the Appendix 2. Fig. 22 shows the results 
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of the calculation for the minimum detectable degree of polarization versus incident 
photon energy for a 3 s detection of a 100% polarized source with the spectrum of 
the Crab (10 hr observation). 

References for Section 2.4. 

Aitbaev, F.B. et al, Proc. 17th ICRC 5, 426 (1981). 

Butler, R.C., et al., Nucl. Instr. and Meth. 221, 41 (1984). 

Costa, E., et al., Astrophys. Space Sci. 100, 165 (1984). 

Daniel, R.R., and Stephens, S.A., Rev. Geophys. Space Phys. 12, 233 (1974). 
Dunphy, P.P. et al., Nucl. Instr. and Meth. A274, 362 (1989). 

Owens, A. et al., Proc. 19th ICRC 3, 314 (1985). 

Schonfelder, V. et al, Astrophys. J. 240, 350 (1980). 

Thompson, D., JGR 76, 1309 (1974). 

Yerma, S.D., JGR 72, 915 (1967). 


2.5 Design of 10 Liter LXe-TPC/Electronics/Gas System 

The only convincing evidence of the suitability of the LXe-TPC technology for 
gamma-ray observations in space is to demonstrate the technology with a laboratory 
prototype of similar size as that of the proposed flight instrument. With substantial 
contribution from the Waseda collaborators, the design of a 10 liters LXe-TPC 
prototype was initiated during 1993. In its final version, the LXe detector and its 
cryostat is shown in Fig. 23. The internal electrode structure of the TPC is shown 
schematically in Fig. 24. For cost reasons, the sensitive area and depth of liquid is 
half the scale of the future balloon version. Fig. 25 shows the gas purification and 
recovery system associated with the 10 liter detector. Its construction and test were 
completed in 1993. Xe gas equivalent to five liquid liters were already procured with 
1993 funds and purified with such system. The Xe is stored in two SS cylinders 
(each of 4 gallons capacity) . The design of the analog and digital readout electronics 
for the 10 liter prototype was also done during 1993. A schematic of the readout 
chain is shown in Fig. 26. The induction signals (50 from X-wires and 50 from 
Y-wires) and the collection signal from the anode are directly coupled to low noise 
charge sensitive preamplifiers. Additional signal amplification needed to drive the 
FADCs is obtained with a second stage amplification system which also include a 
band pass filter for noise reduction. The analog signals are then digitized with a 
VME based FADC system. Both the amplitute and the timing of each event is 
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recorded. For more details on the characteristics of the electronics components we 
refer to the Proposal Renewal for 1994, submitted separately to NASA. 
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A Monte Carlo analysis of the liquid xenon TPC as gamma-ray 
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E. Aprile, A. Bolotnikov, D. Chen and R. Mukherjee 
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Extensive Monte Carlo modeling of a coded aperture y-ray telescope based on a high resolution liquid xenon TPC has been 
performed. Results on efficiency, background reduction capability and source flux sensitivity arc presented. We discuss in 
particular the development of a reconstruction algorithm for events with multiple interaction points. From the energy and spatial 
information, the kinematics of Compton scattering is used to identify and reduce background events, as well as to improve the 
detector response in the few MeV region. Assuming a spatial resolution of 1 mm RMS and an energy resolution of 4.5% FWHM at 
I MeV, the algorithm is capable of reducing the background rate expected at balloon altitude by an order of magnitude, thus 
significantly improving the telescope sensitivity. 


1. Introduction 

Significant progress in experimental y-ray astro- 
physics can only be achieved with the development of 
instruments with improved imaging capability. Efforts 
to improve the sensitivity of a y-ray telescope without a 
corresponding improvement in angular resolution will 
result in source confusion in the ficld-of-vicw (FOV) 
and will limit the mapping of extended source distribu- 
tions. Two important scientific objectives which typify 
the need for detectors with improved imaging capabil- 
ity, along with good sensitivity and energy resolution, 
arc the 1.8 MeV line emission from radioactive : *AI 
and the 0.511 MeV positron -electron annihilation 
emission. For both cases a precise study of the spatial 
distribution of the emission will uniquely identify the 
source of radiation. Among the proposed novel tech- 
niques for imaging astrophysical y-ray sources in the 
low to medium energy region, the liquid xenon time 
projection chamber (LXc-TPC) is recognized as very 
promising and worth a strong R&D effort. Like an 
electronic bubble chamber, a TPC is capable of visual- 
izing the complex histories of y-ray events with multi- 
ple interactions, initiated by either Compton scattering 
or pair production. As a result of this imaging, efficient 
background rejection is also achieved, reducing the 
requirement for a massive anticoincidence shield of the 
type that is required for germanium and sodium iodide 
detectors. From the energy and spatial information 
available for each point of interaction, Compton kine- 
matics enables the reconstruction of the direction of 
the incoming y-ray. The direction is not unique unless 


the Compton scattered electron is also imaged. Given 
the short range of low energy elect rons in liquid xenon, 
this is not possible for a practical spatial resolution of 
the order of millimeters. In the LXc-TPC Compton 
telescope, as proposed by us in ref. fl), the angular 
resolution for the MeV region is limited by the small 
separation between two successive y-ray interactions. 

To overcome this limitation and release the strong 
requirement on the TPC spatial resolution, wc propose 
to combine the imaging capability of the LXc-TPC 
with that of the coded aperture technique, for studying 
astrophysical sources in the few' MeV region with 
arc min accuracy. In a coded aperture telescope, a 
mask consisting of a two-dimensional array of opaque 
and transparent elements is placed between the source 
and a position sensitive y-ray detector. Every source 
within the FOV casts a shadow of part of the mask 
onto the detector. By properly decoding the pattern 
obtained at the detector plane, an image of the source 
is inferred. The successful application of coded aper- 
ture imaging in the low energy y-ray region has been 
demonstrated hy two balloon borne telescopes: the 
New' Hampshire Directional Gamma-Ray Telescope, 
DGT [2], the Caltech Gamma-Ray Imaging Payload, 
GRIP [3], and more recently the SIGMA telescope 
aboard the GRANAT satellite [4], All these telescopes 
use scintillators as position sensitive y-ray detectors. 
The advantage of the LXc-TPC/codcd telescope is its 
unique capability to suppress background over a wide 
energy range, thus improving the sensitivity to weak 
sources, while maintaining an excellent detection effi- 
ciency, energy and angular resolution. 
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In this paper the technical aspects of the LXe-TPC 
development are not considered. For this we refer to a 
previous publication by Aprile et al. [5]. We present 
here initial results from a Monte Carlo analysis of the 
expected performance of the LXe-TPC coupled to a 
coded aperture as y-ray telescope for 0.3-10 MeV. 


2. Monte Carlo analysis and results 

2.7. Telescope model 

The mode of the y-ray telescope used in the Monte 
Carlo simulation is schematically shown in fig. 1. It 
consists of a coded aperture mask, located 1 m above 
the LXe-TPC y-ray detector. An active Csl shield has 
been included. In practice, the background rejection 
capability of the TPC reduces the need for massive 
anticoincidence shield. The sensitive area of the TPC is 
39 x 28 cm 2 . The total depth of liquid xenon is 10 cm. 
The detector's dimensions and mask-detector separa- 
tion have been determined on the basis of the existing 
gondola for the DGT telescope which we plan to use 
for the first balloon flight of the telescope [6]. The TPC 
is operated in the ionization mode. Both the energy 
and the spatial distribution of each ionizing event 
within its sensitive volume is measured. The coordi- 
nates in the X-Y plane arc inferred from the signals 
induced on two orthogonal wire planes, while the coor- 
dinate along the Z-dircction is inferred from the drift 
time measured with respect to a zero time provided by 
the fast scintillation signal. As indicated by our first 


experimental results on the induction signals in a 3.5 1 
LXe-TPC, to be published elsewhere, the spatial reso- 
lution in the plane is better than 5/yT 2, where S is 
the spacing between the induction wires. For the spa- 
tial resolution on the coordinate along the drift axis we 
have previously measured a, = 0.18 mm [5]. For the 
Monte Carlo calculations, a resolution of <r x * <r y — I 
mm and a z — 0.2 mm has been assumed, along with the 
experimental energy resolution of 4.5% FWHM at 1 
MeV, extrapolated according to a E~ ]/2 law {7J. 

The mask is a 2x2 mosaic of a basic uniformly 
redundant array (URA) [8]. It consists of a 43 X 41 
element pattern of 0.91 cm X 0.68 cm X L2 cm thick 
tungsten blocks. With the assumed separation of 1 m 
between mask and the detector, an angular pixel size 
of 0.65° X 0.49° and a FOV of 28° X 20° is defined. 
Much finer angular resolution is obtained by increasing 
the mask-dctcctor separation. The precision in point 
source localization depends, however, also on the de- 
tector’s spatial resolution and the statistical signifi- 
cance of the source data. We have estimated that a 
precision at the areminute level is possible with a 
practical spatial resolution of 1 or 2 mm. 

2.2. Reconstruction algorithm based on Compton kine- 
matics 

The interaction of y-rays in liquid xenon is domi- 
nated by Compton scattering and pair production, for 
energies larger than few hundred kcV. The typical 
cvcnC registered by the TPC is therefore characterized 
by multiple interaction points before the original y-ray 



Fig. 1. Schematic of the LXe-TPC/coded mask y-ray telescope. 
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Fig. 2. Distribution of interaction points for 7 -ray energies of 
(a) 500 keV, (b) 1 MeV and (c) 5 MeV. 


energy is fully deposited in the liquid. This is illus- 
trated in fig. 2 where the distribution of the number of 
interaction points for three different y-ray energies is 
plotted. The probability that the first interaction point 
is a Compton scattering is 90% at 1 MeV, decreasing 


to about 25% at 10 MeV. To improve the signal/ noise 
level in the few MeV region where multiple Compton 
events dominate, an algorithm based on the kinematics 
of the Compton process was developed and tested. The 
algorithm uses the capability of the LXe-TPC to visual- 
ize each incoming y-ray event and to reconstruct the 
initial direction from the energy deposited in each 
interaction point and the coordinates of the points. 
Consider an event in which a y-ray produces N energy 
deposition points. Let W Q be the energy of the incom- 
ing photon, the energy of the scattered photon, E i 
the energy of the scattered electron and cos 0, the 
cosine of the scattering angle. If the coordinates and 
the energy depositions for each interaction point as 
well as the right order of the successive interactions 
are known, then from energy and momentum conserva- 
tion we can write for each interaction point: 


1 - cos Q t ~ l/W t - \/W t _ ,, (1) 

where W i = W t ± 1 + E 4> , , for / = 1 , N - 1 . For any given 
distribution of interaction points generated inside the 
sensitive volume, we now assume that the last interac- 
tion point is a photo-absorption and that the initial 
direction of the incoming y-ray is known. Calling 

1 - COS 0 , a * 


eq. (1) becomes 
W,^W t =a t . 

Similarly, for the (/ + l)th point. 
From cqs. (3) and (4), 


W, = 


e, + e i41 


(3) 

(4) 

(5) 


Also, if the total energy is contained, W {) is 


^0- 


<*i 

E 2 + Ey + ■ ' • E n 


( 6 ) 


Alternatively, from energy conservation alone wc have 
an independent estimate of the energy of the scattering 
photon at each point. Wc express this as IV/. 

H/ = E l+t + E„ (7) 

for 1-0, N - 1. To test the validity of the assumption 
that the total energy is contained, as well as if the 
assumed initial direction of the y-ray is kinematically 
possible, the function 

/V- 1 

r = l w- «:t (8) 

/«<) 

is minimized for all possible combination of points. 
The calculation takes into account the errors on the 
energy and spatial coordinates. 






E Aprile et ai / Monte Carlo analysis of liquid Xe TPC 


219 




Fig. 3. Efficiency of the LXe-TPC for normally incident 
photons. 


Fig. 5. Ratio of events under the peak to that under the 
Compton continuum for the LXe-TPC/coded mask before 
and after event reconstruction. 


2.3. Detection efficiency 

The Monte Carlo calculations were based on the 
EGS4 computer code [9], The excellent 7-ray detection 
efficiency of liquid xenon over a wide energy range is 
shown in fig. 3. The full energy peak efficiency was 
calculated for 7-rays in the range 0.3-10 MeV, nor- 
mally incident on the forward face of a 10 cm deep 
liquid xenon detector. The energy of the events in the 
sample was required to be totally contained. The use of 
a coded mask reduces the detection efficiency by a 
factor of two, assuming 100% opacity of the closed 
mask elements. This is shown in the upper curve of fig. 
4. The effect of the 3 mm stainless steel detector wall 
was also accounted for. The effect of the Compton 
reconstruction algorithm, described in section 2.2, on 
the telescope efficiency is also shown in fig. 4 as a □. 



Fig, 4. Efficiency of the LXe-TPC/coded mask, before and 
after event reconstruction. 


Because the algorithm is tailored for events which start 
with a Compton scattering, all events with a single 
photo-absorption or starting with a pair production are 
rejected. To recover the detection efficiency at low and 
high energy, different reconstruction techniques can be 
applied. For example, the majority of low energy pho- 
tons coming from the source will interact in the top 
region of the LXc-TPC. Fiducial volume cuts for single 
site events can therefore drastically improve the sig- 
nal/ noise ratio. This is shown in fig. 4 as a •. All 
single-site events reconstructed outside the first 1 cm 
layer of liquid xenon from the top were rejected. For 
high energy events which start with a pair production, 
tracking of the electron/ positron can be used on an 
event -by-eve nt basis, to estimate the 7-ray direction 
from the opening angle of the pair. Such an algorithm 
has however not been applied as one can sec from the 
drastic drop in efficiency at energies above ~ 2 McV. 

The capability of the reconstruction algorithm to 
reject events under the Compton continuum is shown 
in fig. 5, which shows the ratio of the total counts 
under the peak to that under the Compton continuum. 
A large improvement in this ratio is achieved for the 
MeV region. The excellent capability of the same algo- 
rithm to recognize 7-ray events entering the TPC nor- 
mally from the bottom was also tested. Results arc 
shown in fig. fi. For this calculation total energy con- 
tainment was not required. 

2.4. Background reduction and flux sensitivity 

The ultimate flux sensitivity is perhaps the most 
relevant “figurc-of-merir of a 7-ray telescope. Among 
the factors which determine the sensitivity are: i) the 
energy of the 7-rays and ii) the detector’s effective 
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Fig. 6. Percentage of events entering the LXe-TPC from the 
bottom before and after event reconstruction. 


area, energy and spatial resolution, the exposure time 
and the total background rate. To calculate the back- 
ground expected in the LXe-TPC/coded mask tele- 
scope at balloon altitudes, we have taken into account 
the dominant atmospheric and cosmic diffuse compo- 
nents entering the forward aperture of the telescope or 
leaking through a thick (5 cm) Csl active shield. The 
flux and angular distribution of the atmospheric y-rays 
used in the calculation were taken from the parameter- 
ized forms given by Costa ct al. [10] and the cosmic 
diffuse spectrum used was that given in ref. [11]. The 
internal backgrounds from natural radioactivity, cosmic 
ray induced radioactivity and activation of instrument 



Fig. 7. Calculated background count rate before and after 
event reconstruction. 

materials have been neglected, as we expect to reject 
the majority of these single site events. Pulse shape 
discrimination of the scintillation light pulses offers 
additional background rejection capability, as it is 
strongly dependent on the ionization density of the 
event. For example, neutron induced events can easily 
be discriminated against, based on the light signature. 
Fig. 7 shows the calculated background flux, before 




Fig. 8. The 3 <r minimum (a) line and (b) continuum flux sensitivity of the LXe-TPC/coded mask 7-ray telescope. 





E. Aprile et al. / Monte Carlo analysis of liquid Xe TPC 


221 


and after applying the Compton event reconstruction 
of section 2.2. and the fiducial volume cut. We as- 
sumed a float altitude of 3 g/cm 2 over Palestine, 
Texas. The integrated flux over the 0.1-10 MeV region 
gives about 340 counts/s, consistent with typical back- 
ground rates measured at these altitudes. A back- 
ground reduction of almost a factor of 10 is obtained 
by identifying 7-rays which kinematically could not 
have come through the FOV of the telescope. It is 
clear that even better background suppression can be 
achieved, on an event-by-event basis, using different 
algorithms for different 7-ray energies, as previously 
discussed. 

From the calculated efficiency and background rate 
we have estimated the 3<r minimum flux sensitivity of 
the proposed telescope. The results are shown in figs. 
8 a and 8b. For a typical exposure of 3 x 10 4 s, we find 
a value of 9 x 10 -5 photons/cm 2 s for 511 keV and of 
6xl0“ 5 photons/cm 2 s for 1.8 MeV. A continuum 
sensitivity of 3 X 10" 7 photons/cm 2 s is obtained at 1 
MeV. The sensitivity curves of the other instruments, 
shown for comparison, have been taken from Winkler 
[12]. When combined with the arcminute source local- 
ization accuracy, the high sensitivity of the LXe-TPC 
telescope makes it competitive with many satellite in- 
struments, even with the much shorter observation 
time available in a balloon flight. 


4. Conclusion 

Monte Carlo modeling of a realistic LXc- 
TPC/coded mask telescope has demonstrated the 
unique capability of this instrument to suppress back- 
ground events as a direct consequence of the imaging 
properties of the LXe-TPC. By using a reconstruction 
algorithm for multiple site events, based on the kine- 
matics of Compton scattering, wc find that the ex- 
pected background flux at balloon altitude is reduced 
by nearly an order of magnitude, resulting in superior 
line flux sensitivity in the MeV region, compared to 


that of existing scintillator-coded mask or proposed 
segmented germanium/ coded mask 7-ray telescopes. 
It is clear that this is only a first step towards the 
application of the liquid xenon technique for precise 
imaging of cosmic 7-ray sources. The full potential will 
be realized when a detector with large effective area, 
improved energy resolution and sub-millimeter spatial 
resolution becomes practically operational as a Comp- 
ton/pair production 7-ray telescope. 
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The Polarization Sensitivity of the Liquid Xenon Imaging Telescope 
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Abstract 

The properties and the expected performance of a liquid xenon (LXe) 7-ray imaging 
telescope, optimized for the MeV energy region are presented. The unique potential of 
this telescope as a Compton polarimeter is particularly emphasized. Based on Monte 
Carlo simulations we show that the modulation factor is as high as 40% at 1 MeV with 
a detection efficiency close to 20%. These figures of merit combined with the excellent 
background suppression capability of the three-dimensional position sensitive LXe detector, 
yield sensitivity at the three sigma level to polarization fractions as small as a few percent 
for strong sources, even in a balloon flight. 

Subject headings: 7 -rays; instruments; polarization 


1. Introduction 

Future missions in low and medium energy 7 -ray astronomy will rely on new telescopes 
with superior imaging capability, high energy resolution, wide sky coverage and high sen- 
sitivity. Of the new detector technologies that have been proposed to meet with these 
requirements, a Liquid Xenon Time Projection Chamber (LXe-TPC) (Aprile, M ukh erjee 
& Suzuki 1989) is among the most promising. The properties of liquid xenon make it an 
excellent radiation detector medium, particularly for 7 -rays. When used in an ionization 
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chamber, operated in the time projection mode, LXe offers an ideal combination of high 
detection efficient, 3-D event imaging with submillimeter spatial resolution and very good 
energy resolution response. Like an electronic bubble chamber, a LXe-TPC is able to visu- 
alize any ionizing event occurring within its sensitive volume. 7-ray events with multiple 
Compton interactions are recognized as such, thus substantially increasing the detection 
efficiency in the difficult Compton region. The interaction of 7-rays in most materials is 
dominated by Compton scattering for energies larger than few hundred keV. As shown in 
Figure 1, the probability that a 7-ray interacts first via Compton scattering in a 10 cm deep 
LXe detector varies between 50% at 1 MeV and 10% at 10 MeV. Also shown is the fraction 
of events with two or more consecutive Compton interactions. In the few MeV region this 
fraction is not small. These multiple Compton events which have a high probability to be 
rejected in a conventional double scatter Compton telescope, carry as much information 
as single Compton events, if all interaction points and the deposited energy in each point 
are measured. The imaging capability can be used not only to analyze the topology of 
the events under investigation, but also to reject background. Since the Compton scatter- 
ing cross-section depends on the polarization of the incident 7-ray, a direct consequence 
of the LXe-TPC intrinsic imaging capability is its sensitivity as a Compton polarimeter. 
To estimate this polarization sensitivity a Monte Carlo program that explicitly considers 
the polarization of Compton scattered 7-rays was written. Results on modulation factor, 
detection efficiency and background rate in the 0.3-10 MeV energy range are combined to 
obtain the polarization sensitivity curve for a 100% polarized source with the spectrum of 
the Crab. 


2. Telescope description and properties 


Figure 2 shows the schematic of the LXe-TPC 7-ray imaging telescope which we have 
proposed for the observation of a variety of discrete and diffuse 7-ray sources in the MeV 
region (Aprile et a 1 . 1992a; Aprile et 3.1. 1993a). It consists of a 3-D position sensitive 
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LXe 7-ray detector, an active veto shield, and a coded aperture mask. The coded mask 
consists of a 2 x 2 mosaic of a basic Uniformly Redundant Array (URA), with a 85 x 83 
element pattern of 0.91 x 0.58 x 1.2 cm 3 thick blocks of tungsten alloy. With 1 meter 
distance between the mask and the detector plane, the FWHM angular resolution is ~ 30 7 
and the nominal fully coded FOV is 28° x 20° FWHM. For a Hkr source strength, the point 
source localization accuracy is estimated to be ~ 1 arcminute, based on the LXe detector 
spatial resolution of 1 mm RMS. The LXe detector has an active area of 39 x 28 cm 2 , and 
a thickness of 30 g cm -2 which gives good sensitivity up to 10 MeV. The lowest energy 
threshold will be determined mostly by the S/N ratio achievable in a practical detector 
operated in the ionization mode. The LXe-TPC works on the principle that the free 
ionization electrons liberated in the liquid by a charged particle can drift, under a uniform 
electric field, towards a signal readout structure. The ionization signals induced on the 
sensing elements provide both the spatial and total energy information for each event. In 
our design the information in the XY plane is obtained from the signals induced on two 
orthogonal wire planes while the Z information, along the direction of drift, is inferred 
from the drift time, measured with respect to a zero time provided by the fast signal from 
the primary scintillation in LXe. The technical feasibility of such a detector has been 
studied at Columbia for the past few years (Aprile et al. 1992a, 1992b; Aprile, Mukherjee 
& Suzuki 1991a, 1991b, 1990). We have finished the construction and are currently testing 
with laboratory 7-ray sources a 10 liters detector, with approximately half the sensitive 
area and half the drift space as the one proposed for the first balloon flight. 

3. Detector Performance: Monte Carlo Results 


The Compton scattering process is an unique tool for determining the polarization of 
the incident 7-ray because there is a large difference in the scattering cross-section for the 
orthogonal states of polarization. The relation between the energy of the incident 7-ray, 
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Eq, and that of the scattered 7 -ray, E, is given by: 


I? = £b [l + -Eo (l — cos#)] 1 


( 1 ) 


where E and Eq are expressed in units of moc 2 , and 0 is the scattering angle. If the 
incident 7-ray is completely polarized, then the differential cross section for being Compton 
scattered is: 

** “ \ r «' I []| + % - 2sin 2 tfcos 2 ^]<« (2) 

where r 0 is the classical radius of the electron, <j> is the azimuthal angle, and dfl is the 
element of solid angle into which the 7 -ray is scattered. The polarization sensitive part 
of the cross section is contained in the cos 2 4> term. A measure of the response of the 
Compton scattering process to polarized radiation is obtained from the “asymmetry ratio”, 
R, defined as 

R = dagojdcro ( 3 ) 


where da 0 and dago are the differential scattering cross sections for 0 = 0° and <j> = 90°. 
Figure 3 shows the asymmetry ratio as a function of 6 for different incident 7 -ray energies. 
At low energies, for values of 6 close to 90°, the Compton scattering process provides ideal 
response to radiation polarized in the reaction plane. As the energy of the incident 7 -ray 
increases, the angle 6 at which R is maximum, # max , decreases. 

The modulation factor of a Compton polarimeter is defined as: 


W) = 


JVx(#) + JV„(#) 


(4) 


where N±(6) is the detected count rate of the scattered 7 -ray in a direction perpendicular 
to the 7 -ray’s electric vector (that is, for <j> = 90°) and N\\ (6) is the detected count rate 
in the direction parallel to the electric vector (that is, for <f> = 0°). Q is a measure of the 
suitability of the detector as a Compton polarimeter. The advantage of the LXe-TPC over 
more conventional double scatter polarimeters is that the scatterer and absorber are the 
same detector, which largely increases the detection efficiency, with Aix acceptance. 
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To evaluate the performance of the LXe telescope as a Compton polarimeter we have 
written a Monte Carlo program based on EGS4, which explicitly treats the polarization 
dependence of the Compton process. For the simulations we assumed a LXe active volume 
of 39 x 28 x 10 cm 3 and a spatial and energy resolution of: a x = <j y = 1 mm, a z = 0.2 mm 
A E/E = 4.5% FWHM at 1 MeV with a E~ o s dependence. A flux of 100 % polarized 7 -rays 
incident normally on the detector plane was generated. As previously mentioned, a large 
fraction of these events are characterized by multiple Compton scatterings. For each event 
the detector measures the coordinates and the energy deposited in each interaction point. 
FVom the locations of the first two successive interaction points the scattering angle 6 and 
the azimuthal angle <f> of the scattered 7 -ray is inferred, provided that the correct order of 
interaction points is known. If the energy is totally contained, the original 7 -ray energy 
is also inferred. To find the correct order of the interactions, a reconstruction algorithm 
based on Compton kinematics was developed (Aprile et a 1. 1993b). The efficiency of 
the LXe polarimeter lies in its capability to include multiple Compton events which in a 
conventional double scatter polarimeter are rejected. 

As shown in Fig. 2 , the detector is most sensitive as a polarimeter for a particular 
scattering angle, 0 max , when the asymmetry ratio is a maximum. To increase the detector 
count rate, all events were accepted for which the scattered 7 -ray was within some angle 
bin about 6 m&x . Figure 4 shows the distribution of the scattered 7 -rays in the range <f> = 0 ° 
to <f> 90 , for different A 0 bins, where A 6 = \{0 $rnax)|- The simulations were performed 

for 100% polarized 7 -rays of energy 500 keV, incident normally on the detector surface. 
The results for the modulation factor, Q, versus incident 7 -ray energy are shown in Fig. 

5, for A0 - 5°. At 1 MeV Q is as high as 40%. For comparison, the modulation factors 
for the Ge(Li) polarimeter (Taras 1968), the Si(Li)-Ge polarimeter (Ohya et al. 1989), 

V 

and the CsI(Tl) polarimeter (1 cm 2 cross-sectional area) of the Imager on INTEGRAL 
(Swinyard et al. 1991) are also shown. 

The polarization sensitivity of a polarimeter, or the minimum detectable degree of 
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polarization, for a given source flux and observation time, is determined by the modulation 
factor, as well as by the detection efficiency and the background rate. The polarization 
sensitivity, P min of the LXe-TPC polarimeter was calculated at the 3 a level for a 100% 
polarized source with the spectrum of the Crab, by using the relation: 

p3<r ^ / c/gpuroe 4~ B 

mlD ^source V AEAt ^ 

where € is the detector efficiency, Q is the modulation factor, J TOUrce is the source flux, B 
is the background rate, t is the observation time, A is the sensitive area of the detector, 
and A E is the energy bin. 

The Crab spectrum for the MeV region was taken from the DGT experimental results 
as 5.1 x 10- 3 P- 188 cm -2 s“ 1 MeV~ 1 (Dunphy et al. 1989). The efficiency, and background 
rate at balloon altitude, were estimated using Monte Carlo simulation (Aprile et al. 1993b) 
and are shown in Fig. 6 and Fig. 7. We note that although the use of the coded mask 
reduces the detection efficiency by a factor two, a 20% efficiency is still achieved at 1 MeV. 
For the background estimation, we considered only the dominant atmospheric and cosmic 
diffuse components and assumed an active 5 cm thick Csl anti-coincidence shield. As 
shown in Fig. 7, a background reduction of almost a factor often is achieved by correctly 
identifying 7-rays which kinematically could not have come through the telescope’s FOV. 
The background rate, after event reconstruction, as used in the polarization sensitivity 
calculation, was approximated by a power law as 10 x P _2 s _1 MeV -1 . The energy bin A E 
was set equal to 200 keV and the observation time was taken as 10 h, for a typical balloon 
flight. The results of the calculation are shown in Fig. 8. The histogram stops at 2 MeV 
since, at higher energies the polarization of the Crab source cannot be measured at the 
3 a level, even if we assume a 100% polarized flux. 
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4. Conclusions 


Based on Monte Carlo simulations we conclude that the LXe-TPC imaging telescope 
will work as a sensitive Compton polarimeter for cosmic 7-ray sources in the MeV region, 
even as a balloon payload. A modulation factor of 40% was calculated at 1 MeV. This 
combined with a detection efficiency close to 20% and the low background rate which 
can be achieved with an intrinsically imaging 7-ray detector, will permit sensitivity to 
polarization fractions as small as few percent for strong sources. 

This work was supported by NASA grant NAGW-2013. 
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Figure Captions 

Fig. 1 . Expected number of 7-ray events with a single Compton scattering and with two or 
more consecutive Compton scatterings, vs. incident 7-ray energy. 

Fig* 2. Schematic of the LXe-TPC/coded mask imaging 7-ray telescope. 

Fig. 3 . Asymmetry ratio, R, as a function of the scattering angle 6 . 

Fig. 4 . Distribution of the scattered 7-rays for <f> = 0° - 90 °, simulated for 100% polarized 
incident 7-rays of energy 500 keV. 

Fig. 5 . Modulation factor, Q, vs. incident 7-ray energy for the LXe-TPC. For comparison, 
results on Q for other polarimeters are also shown. 

Fig. 6. Detection efficiency vs. incident 7-ray energy, before and after event reconstruction. 
Fig. 7 . Background rate (atmospheric and cosmic diffuse components only) before and after 
event reconstruction. 

Fig. 8. The minimum detectable degree of polarization vs. incident photon energy for a 3 <7 

detection of a 100% polarized 7-ray source with the spectrum of the Crab for an 
observation time of 10 hours. 
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A Liquid Xenon Imaging Telescope for Gamma-Ray 
Astrophysics: Design and Expected Performance 

E. Aprile, A. Bolotnikov, D. Chen and R. Mukherjee 
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A high resolution telescope for imaging cosmic 7-ray sources in the MeV region, with an angular 
resolution better than 0.5° is being developed as balloon-borne payload. The instrument consists of a 
3-D liquid xenon TPC as 7-ray detector, coupled with a coded aperture at a distance of 1 meter. A study 
of the actual source distribution of the 1.809 MeV line from the decay of 26 Al and the 511 keV positron- 
electron annihiliation line is among the scientific objectives, along with a search for new 7-ray sources. 
The telescope design parameters and expected minimum flux sensitivity to line and continuum radiation 
are presented. The unique capability of the LXe-TPC as a Compton Polarimeter is also discussed. 


1. INTRODUCTION 

Gamma-ray telescopes with true imaging ca- 
pability and high flux sensitivity are essential for 
studying the highest-energy phenomena in the 
universe. Fine imaging provides accurate posi- 
tioning of the sources detected within the FOV 
and good angular resolution to map regions of 
diffuse emission and separate point source contri- 
butions. The importance of true source imaging 
is particularly evident in the study of two of the 
most pressing problems in low energy 7-ray as- 
tronomy: the 1.089 MeV line emission from the 
decay of 26 A1 and the 511 keV positron-electron 
annihilation line emission from the Galactic Cen- 
ter. 

In 1977, Ramaty and Lingenfelter [1] sug- 
gested that galactic nucleo-synthetic production 
of 26 AI in supernova events over the past few mil- 
lion years could give rise to a detectable 7-ray 
line at 1.809 MeV. This line arises from the elec- 
tron capture (18%) or positron decay (82%) of the 
million-year mean life 26 A1 and was first detected 
in 1984 [2] at a flux level of 4.3±0.8 xlO -4 pho 
tons cm' 2 s' 1 rad' 1 at the Galactic Center. Sev- 
eral subsequent confirmations of the line energy 
and flux level have been made. Some potential 
sources of 26 A 1, which have been proposed, are 
supernovae, novae, red giants in the Asymptotic 
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Giant Branch (AGB), Wolf-Rayet stars or nearby 
OB stars (see e.g., [3] for a recent review). Since 
these objects have more or less known or inferred 
galactic distributions, it is believed that a mea- 
surement of the spatial distribution of the 26 AI 
1.809 MeV line intensity will identify the 26 A1 
source. The only instrument which could mea- 
sure this radiation with imaging capability, is the 
Compton telescope, the most advanced version of 
which is COMPTEL on the COMPTON Obser- 
vatory. COMPTEL however cannot directly mea- 
sure the 1.809 MeV spatial distribution. The only 
definite statement that can be made about the 
26 A1 spatial distribution from the latest COMP- 
TEL results, at the present time, is that a point 
source near the Galactic Center can be excluded 
[4]. Clearly, there is a requirement to measure 
directly the spatial distribution of the 1.809 MeV 
line with a true imaging telescope. 

As for the 511 keV line, the debate between 
point like and diffuse nature of the emission con- 
tinues to date and can only be fully resolved with 
a high level imaging map of the Galactic Center 
region at 7-ray energies. 

At the present time the OSSE instrument 
on the COMPTON Observatory is mapping the 
distribution of the annihilation line [5]. Since 
the OSSE measurements give the lowest galac- 
tic center flux measurements so far Skibo, Ra- 
maty and Leventhal [6] have used these results 
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an d other off-center measurements to test differ- 
ent models for the origin of the diffuse or steady 
galactic plane 511 keV component- On the other 
hand, the origin of the variable narrow line galac- 
tic 511 keV radiation may be associated with the 
bright hard X-ray source 1E1740.7-2942 which 
was studied by the imaging telescope SIGMA 
on the GRAN AT satellite during the spring-fall 
of 1990 and in early 1991. Sunyaev tt ai [7] 
have identified three spectral states for this source 
which range from a “low state” a normal (Cygnus 
X-l like) state to a hard state in which a bump 
appears in the spectrum between (300-600) keV. 
The broad feature of the spectrum has been in- 
terpreted as annihilation of positrons in a hot 
medium (~ 40 keV). This is consistent with the 
temperature of the accretion disk derived from 
the X-ray continuum spectrum. 

Subsequently it was proposed [8,9] that in 
addition this high energy source injects positrons 
into a molecular cloud where they slow down and 
annihilate to produce the narrow component of 
the 511 keV line emission. 

Future studies of the 511 keV emission re- 
quire the most advanced imaging telescope with 
good to excellent energy resolution. 

Of the techniques proposed for 7-ray imag- 
ing and spectroscopy of astrophysical sources, the 
Liquid Xenon Time Projection Chamber (LXe- 
TPC) is among the most promising. The prop- 
erties of liquid xenon make it very efficient for 7- 
ray detection. When used in an ionization cham- 
ber, operated in the time projection mode, this 
medium offers a combination of high detection 
efficiency, excellent spatial resolution and very 
good energy resolution. Like an electronic bub- 
ble chamber, a LXe-TPC with three-dimensional 
position sensitivity is capable of visualizing the 
complex histories of 7-ray events initiated by ei- 
ther Compton scattering or pair-production. As 
a result, efficient background rejection is also 
achieved, reducing the requirement for massive 
anticoincidence shielding of the type that is re- 
quired for germanium or sodium iodide 7-ray de- 
tectors. The angular resolution of the LXe-TPC 
as a Compton telescope is however limited, in the 
few MoV region, by the small separation between 
two successive 7-ray interactions [13], To achieve 


imaging with good angular resolution at low en- 
ergies, the combination of the imaging LXe-TPC 
with a coded aperture is proposed. 

A unique consequence of the LXe-TPC imag- 
ing capability is its sensitivity as a Compton po- 
larimeter. Besides the precise determination of 
the energy and incident direction of a photon, 
determination of its polarization state can give 
further information on the source of 7 rays. The 
main production mechanisms which can give po- 
larized 7 rays are: bremsstrahlung from elec- 
tron beams, electron synchrotron radiation, elec- 
tron curvature radiation, and 7 rays from de- 
excitation of nuclei excited by directed ion beams. 
In the case of the Crab Nebula it has been deter- 
mined that the nebular X-ray emission is polar- 
ized [10]. Existence of UHE (> 10 14 ) electrons in 
this source could yield polarized nebular 7-rays of 
a few MeV. If curvature radiation from electrons 
is the source of MeV 7 rays in pulsars, such as 
the Crab and Vela, then polarization might also 
be expected. 

In general it has been recognized in the study 
of X-ray sources, that measurement of the direc- 
tion and magnitude of the photon polarization 
could significantly contribute to a better under- 
standing of the physical processes in compact ob- 
jects, such as pulsars, Black Holes and AGN. 

2. TELESCOPE DESIGN 

2.1. Introduction 

The telescope is schematically shown in Fig. 
1. It consists of a coded aperture mask, located 
1 meter above a LXe-TPC. The sensitive area of 
the TPC is 39 x 28 cm 2 . The active depth of liq- 
uid xenon is 10 cm. Fig. 2 shows the LXe-TPC 
in more detail. The event trigger to the readout 
electronics is provided by the fast primary scintil- 
lation light detected by two UV sensitive PMTs. 
The intrinsic instrumental angular resolution in 
the coded mask configuration is determined by 
the size of the mask unit cell, the mask-detector 
separation, and by the accuracy to which the pho- 
ton interaction points in the detection plane can 
be determined. The coded mask that we have 
assumed in our design and Monte Carlo simu- 
lations consists of a 85 x 83 element pattern of 
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2 s 1 


Min Flux (Continuum) 3 x 10” 7 ph cm -2 s~ 1 keV 


3a at 1 MeV 


(3 x 10<s) 


Fig. 1. Schematic of the LXe-TPC/coded mask 
imaging 7-ray telescope. 


0.91 x 0.58 x 1.2 cm 3 thick blocks of tungsten al 
l °y* The 1 meter separation between the mask 
and the LXe detector plane is determined by 
the University of New Hampshire gondola [11] 
which we plan to use for the first balloon flight. 
Phis defines a pixel element of angular dimension 
3.65° x 0.49°. The nominal FOV is 28° x 20° 



Fig. 2. Schematic of the LXe-TPC detector. 


and the source localization accuracy is estimated 
to be ~ I areminute, for a 10<7 source strength. 
For background suppression at balloon altitude 
an active shield has been assumed around the de- 
tector. The type and amount of shield needed 
will ultimately be determined by the type of event 
triggering and selection on board, by dead time 
consideration, telemetry rate as well as cost and 
weight consideration . 

2.2. The LXe-TPC: Status of Development 
A LXe-TPC works on the principle that free 
ionization electrons liberated by a charged parti- 
cle in the liquid can drift, under a uniform electric 
field, from their point of creation towards a signal 
read-out region. Here the charge signals induced 
or collected on sensing electrodes are detected to 
yield both the spatial distribution of the ionizing 
event and its energy. 

For 7-rays it is the electrons or positrons cre- 
ated by photoabsorption, Compton scattering or 
pair production, which will ionize as well as ex- 
cite the xenon atoms creating a large number of 
electron-ion pairs and scintillation photons. For 
3-D imaging of 7-ray events in LXe we plan to use 
a sensing electrodes geometry based on the orig- 
inal design by Gatti ei ai [12]. Two orthogonal 
induction wire planes separated from the drift re- 
gion by a screening grid, give the X-Y event infor- 
mation. The measured drift time, referred to the 
scintillation trigger, and the known drift velocity 
provides the Z-in formation. The total event en- 
ergy is measured from the total charge collected 
on an anode plate, placed below the induction 
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wires. 

In order to verify the feasibility of such a 
detector t the Columbia group started in 1989 an 
intensive R&D program on LXe. The attenuation 
length of electrons and UV photons in purified liq- 
uid xenon, the ionization and scintillations yields 
of electrons and alpha particles, the energy and 
spatial resolution have been studied. 

The experimental results obtained on these 
aspects relevant for the development of a liquid 
ionization TPC, are documented in several refer- 
ences [13-18]. Especially relevant are the latest 
experimental results obtained with a 3.5 liter 2D- 
TPC prototype [19] equipped with a multi-wire 
structure to detect the induction signals in liquid 
xenon. The results demonstrate both the capa- 
bility of a large volume LXe detector to provide 
similar or better energy resolution than the pre- 
viously reported value of 4.5% FWHM for 1 MeV 
radiation, as well as the imaging capability. 

Figure 3 $hows an example of collection and 
induction signals produced by a 7-ray event in 
the LXe-TPC prototype. The induction signal, 
which has the expected triangular shape, has a 
large S/N ratio of 12:1, even for a typical point- 
like charge deposition produced by a 7-ray inter- 
action. The dependence of the induced signal on 
the lateral position of the drifting electron cloud 



Fig. 3. Collection signal (upper trace Gain=l) 
and induction signal (lower trace Gain=400) pro- 
duced by a 7-ray event in the 3.5 liter LXe-TPC 
prototype. 


with respect to the wire cell [19], offers the pos- 
sibility to derive the spatial coordinate of each 
event by weighting the signal amplitude on neigh- 
bouring wires. Thus the spatial resolution in the 
X-Y plane, can be better than s/ \Zf2, where s is 
the wire spacing. 

Experimental work on the operation and per- 
formance of the LXe TPC prototype implemented 
for full 3-D imaging and triggered by the scintil- 
lation light is in progress. 

3. TELESCOPE PERFORMANCE: 
MONTE CARLO RESULTS 

3.1. Background Rate and Minimum Flux 
Sensitivity. 

To calculate the background expected in the 
LXe-TPC/coded mask telescope at balloon alti- 
tudes, we have taken into account the dominant 
atmospheric and cosmic diffuse components, en- 
tering the forward aperture of the telescope or 
leaking through the active shield (5 cm thick Csl). 
The flux and angular distribution of the atmo- 
spheric 7-rays used in the calculation were taken 
from the parameterized forms given by Costa ci 
al. [20] and the cosmic diffuse spectrum used was 
that given by Shoiifelder, Graser and Daugherty 
[21]. The internal backgrounds from natural ra- 
dioactivity, cosmic ray induced radioactivity and 
activation of instrument materials have been ne- 
glected, as the majority of these single site events 
can be rejected by simple fiducial volume cuts. , 

The results of the calculation are shown In 
Fig. 4. The integrated flux over the 0.1-10 
MeV region gives about 340 counts/sec, consis- 
tent with typical background rates measured at 
the assumed altitude. An event reconstruction 
algorithm based on the kinematics of Compton 
scattering was developed and used for identifica- 
tion and rejection of background events [22]. As 
shown in Fig. 4, a background reduction of ap- 
proximately a factor of 10 is obtained by iden- 
tifying 7-rays which kinematically couldn’t have 
come through the FOV of the telescope, and by 
applying a fiducial volume cut to remove low en- 
ergy events. 

Based on the calculated 7-ray detection effi- 
ciency [22] and the calculated background rate, 
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Fig. 4. Monte Carlo calculation of the back- 
ground flux at balloon altitude. 


we have obtained the 3c t minimum flux sensi- 
tivity shown in Fig. 5. With a typical balloon 
flight exposure of 3 x 10 4 s, the 3a line sensi- 
tivity is 6 x 10 ~ 5 photons cm -2 s " 1 (1.8 MeV 
line) and 9 x 10 “ 5 photons cm ” 2 s ” 1 (511 keV 
line). The continuum sensitivity is 3 x 10~ 7 pho- 
tons cm ’ 2 s ’ 1 keV ” 1 at 1 MeV. The sensitivity 
curves of the instruments, shown for comparison, 
have been taken from Winkler [23). When com- 
bined with the excellent source localization accu- 
racy, the high sensitivity of the LXe-TPC tele- 
scope makes it competitive with many satellite 
instruments, even with the much shorter obser- 
vation time available in a balloon flight. 

3.2. Simulated Observations of the Crab 
and 511 keV Line 

The Crab Nebula/Pulsar will be the primary 
target for the first verification balloon flight of 
the LXe-TPC coded mask 7 -ray telescope. This 
source is one of the most intense in our energy 



Fig. 5. 3 a minimum line flux sensitivity of the 
LXe-TPC/ coded mask 7 -ray telescope. 

range and is stable, both in intensity and spec- 
trum. Monte Carlo simulations of the expected 
Crab signal were performed using the complete 
telescope system shown in Fig. 1 . 

The Crab Nebula was assumed to be a 
point source in the sky with a spectrum equal 
to 5.5xl0 ” 4 (E/100 keV )” 22 photons cm ” 2 s -1 
keV ” 1 at 10 MeV [24]. The source was aligned 
with the telescope axis and the observing time 
was 10 4 s. The estimated background of Fig. 
4, after event reconstruction, was uniformly dis- 
tributed in the detector’s plane and added to the 
shadowgram of the source. Figure 6 shows the 
resulting deconvolved image of the Crab, for the 
energy interval 0.3 - 0.5 MeV. The Crab signal 
dominates over the background up to several MeV 
with a S/N of about 20a. 

Simulated observations have also been per- 
formed for the low and high state of the 511 keV 
Galactic Center annihilation line. A 10 4 s expo- 
sure time was assumed. The source was placed 
in the center of the FOV, and superimposed on 
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6 . Monte Carlo simulation of the Crab as 
point 7 -ray source in the energy range ( 0 . 3-0 5 
MeV). 



*g- 7. Monte Carlo simulation of the 511 keV 
Galactic Center point source observed for High 
Slate. 


uniformly distributed background of 4xl0 “ 2 
<’« »unts s _1 keV” 1 , as from our estimate. The in- 
tensity of the 51 1 keV line source was chosen to be 
• >x photons cm ” 2 s ” 1 for the “low state” and 
1 * 10“ 3 photons cm ~ 2 s _1 for the “high state”. 


Figure 7 shows the result of the 511 keV image 
for the “high state”. Even in the “low state”, the 
511 keV flux can be detected by our instrument 
at a satisfactory significance level of — 4a. 

3.3. Polarization Sensitivity 

The LXe-TPC imaging capability is also 
ideal to measure the linear polarization of the 
incident 7 -ray undergoing Compton scattering. 
The linear polarization of 7 -rays can be measured 
based on the principle that the Compton scat- 
tering process is sensitive to the polarization of 
the incident 7 -ray, the cross-section for Compton 
scattering being the largest for the case when the 
direction of the scattered 7 -ray is normal to the 
polarization vector of Lhe incident 7 -ray. The ad- 
vantage of a LXe-TPC Compton Polarimeter over 
the conventional Nal(TI), CsI(Tl) or Ge(Li) dou- 
ble scatter Compton telescopes is the enhanced 
detection efficiency offered by a single detector 
working both as scatterer and absorber, as well 
as its combination of good energy and position 
sensitivity. 

A Monte Carlo program was developed to 
estimate the polarization sensitivity of the LXe- 
TPC for a 100 % polarized 7 -ray beam of energy 
varying from 300 keV to 4 MeV, incident nor- 
mally on the detector surface. Figure 8 shows 
the result. For comparison, the polarization sen- 
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Fig. 8 . Monte Carlo calculation of the LXe-TPC 
polarization sensitivity. 
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sitivity of the Ge(Li) polarimeter [25], the Si(Li) 
[26] and the CsI(Tl) polarimeter of the Imager on 
INTEGRAL [27] are also shown. The unique fea- 
ture of the LXe-TPC is its capability to infer the 
scattering angle 0 and the azimuthal angle <j>, with 
an accuracy of about 0.5° [13], for each scattered 
y-ray. We can thus obtain the azimuthal angular 
distribution of the scattered y-rays by selecting 
events from different intervals of scattering angle. 
By applying the detector’s response function, cal- 
culated or measured during calibration tests with 
polarized beams, we can deconvolute the original 
y-ray polarization. Figure 9 shows the modula- 
tion curve in the range <p = 0° to <f> = 90°, sim- 
ulated for 100% polarized y-rays of energy 500 



0 20 <0 60 60 


AN<XC 4 (d*q) 

Fig. 9. Modulation curve of the LXe-TPC for 
100% polarized y-rays of energy 500 keV. 

4. CONCLUSION 

The design and expected performance of a 
y-ray imaging telescope tailored to the 0.3-10 
MeV energy region have been discussed. The 
telescope combines the excellent properties of a 
liquid xenon TPC as 3-D position sensitive y-ray 
detector with the well established imaging prop- 
erties of a coded aperture mask, to achieve high 
efficiency, good spectroscopy and angular resolu- 
tion over the entire energy range of interest. The 
high sensitivity to MeV y-ray lines and contin- 
uum complemented with the good imaging capa- 


bility will permit the observation of a variety of 
astrophysics] sources. Important contributions to 
the field of low energy astrophysics as well as new 
discoveries are expected even in the maiden bal- 
loon flight which is planned for the end of 1994. 

This work was supported by NASA (Award 
NAGW 2013). 
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